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Abstract 

We discuss axisymmetric force-free pulsar magnetospheres with mag- 
netically collimated jets and a disk-wind obtained by numerical solution of 
the pulsar equation. This solution represents an alternative to the quasi- 
spherical wind solutions where a major part of the current flow is in a 
current sheet which is unstable to magnetic field annihilation. 
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1 Introduction 

Interest in the structure of pulsar magnetospheres has been stimulated by Chan- 
dra and Hubble Space Telescope observations of the Crab synchrotron nebula 
which point to an axial-jet/equatorial-disk structure (Hester et al. 2002). An 
analogous structure is observed in the nebula of the Vela pulsar (Pavlov et 
al. 2001). A theoretical model of an aligned rotating pulsar with collimated 
jets was put forward by Sulkanen and Lovelace (1990; SL) who solved the pul- 
sar equation of Scharlemann and Wagoner (1973) on a grid numerically. This 
work was criticized by Contopoulos, Kazanas, and Fendt (1999; CKF) who 
presented numerical calculations showing that the possibly unique solution is 
a quasi-spherical wind without jets but with an equatorial current sheet. The 
quasi-spherical wind solution has been found in the time-dependent, relativistic, 
force-free simulations by Spitkovsky (2004), Komissarov (2006), and McKinncy 
(2006) and in high resolution grid calculations by Timokhin (2006). However, 
the wind solutions may be short-lived owing to the fact that a major part of 
the current flow is in a current sheet which is unstable to magnetic field annihi- 
lation. MHD simulations by Komissarov and Lyubarsky (2004, 2006) indicate 
that a jet-torus configuration can be generated due to the anisotropic energy 
flux density of the pulsar far outside the light cylinder. 
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We reconsider the possibility of jet/disk-wind structures of aligned pulsar 
magnetospheres on the scale of the light-cylinder distance using a different ap- 
proach to the numerical solution of the pulsar equation. We utilize the fact 
that the poloidal current flow along the poloidal field lines within the star's 
light-cylinder is an adjustable parameter. We find that when this parameter 
is sufficiently large, magnetically collimated (±z) jets form within the light- 
cylinder and a quasi-collimated flows outside. The collimation is due to the 
toroidal magnetic field. The analysis by SL did not include the current flows 
outside the light-cylinder and this resulted in a kink in the field lines which 
cross the light-cylinder. In addition to the collimated flows, we find an "anti- 
collimated" disk-wind. The anti-collimation is due to the toroidal magnetic 
field. These solutions have no net poloidal current flow and no current sheets 
inside or outside the light-cylinder. Thus these solutions are not unstable to 
field annihilation. The formation of magnetically collimated jets along the axes 
and an equatorial disk-wind is similar to what is found in the nonrelativistic 
limit for magnetic loops threading an accretion disk (Ustyugova et al. 2000). 
This jet/disk-wind geometry was discussed for the case of pulsars by Romanova, 
Chulsky, and Lovelace (2005). 

Section 2 of the paper discusses the theory, the boundary conditions, and the 
regularity condition at the light-cylinder. It goes on to discuss the conditions 
for having no jets and having jets. Section 3 discusses the numerical solutions. 
For the case of jets we discuss the radial force balance across the jet and the 
vertical force balance in the disk. Section 4 gives the conclusions of this work. 

2 Theory 

The main equations for the plasma follow from the continuity equation V-(pv) = 
0, Ampere's law, V x B = 47rJ/c, Coulomb's law V • E = 47rp e , with p e the 
charge density, Faraday's law, VxE = 0, perfect conductivity, E+v x B/c = 0, 
with v the plasma flow velocity, and the "force-free" condition in the Euler 
equation, p e E + JxB/c = 0. The perfect conductivity implies that E 2 < B 2 . 
Owing to the assumed axisymmctry, — 0, so that the poloidal velocity 
v p = kBj,. Mass conservation then gives B- V(pn) = 0, which implies that pn = 
F(^/)/An, where F is an arbitrary function of the flux function ^ . In cylindrical 
coordinates, B r = —r~ 1 (d^/dz) and B z = r^ 1 (d^ /dr). In a similar way one 
finds that - kB^ = rG(^), so that E = -G(*)V*, and rB^ = #(*), so 
that there are two additional functions, G and H. 

The function G is determined along all of the field lines which go through 
the star. This follows from the perfect conductivity condition at the surface of 
the star, Eg + (v x B)e/c = 0, in terms of spherical (i?, #,</>) coordinates. This 
gives Ee — — (v^Br ~ vrB^/c = —v^Br/c, where vr is zero inside the star. 
Here, = Sl*R*sin9 is the velocity at the star's surface, £1* is the angular 
velocity of the star, and i?* is the star's radius. Thus we have Eg{R st ,9) — 
-Q* (d^/d6)/(R*c), so that G(*) = fi*/c. 

The component of the Euler equation in the direction of V^P gives the force- 



2 



free Grad-Shafranov (GS) equation in cylindrical (r, <f), z) coordinates, 



c 



A** 



c 2 dr 



(i) 



where 



and 



dH{^) 



A* = d 2 /dr 2 - (l/r)(d/dr) + 8 2 /dz 2 , 

(Scharlemann & Wagoner 1973). 

Note that the poloidal current density is given as J p = (c/4 : -jT)'Bp(dH/d i ^) 
where the p— subscript indicates the poloidal component. We consider solutions 
with symmetry about the equatorial plane with for example \&(r, z) — ^(r, —z) 
and B^r, z) = -B<p(r, -z). 

This equation for '5 involves the unknown function J 7 ^) or H(^f). It is 
in general nonlinear. Ampere's law gives § dl ■ B = (Att/c) J dS ■ J, so that 
rB^{r,z) = H(^>) is (2/c) times the current flowing through a circular area of 
radius r (with normal z) labeled by \I/(r, z)= const. 

In the following distances are measured in units of the light cylinder radius, 
tl = c/fi*. The flux function \f is measured in units of \ijr^, where /j, is 
the magnetic moment of the star. The magnetic field is measured in units of 
Bo = mM- 



2.1 Boundary Conditions 

Numerical solutions of equation (1) for &(r, z) are calculated on a uniform grid 
in a region r = to r max and z = to z max with r max = z max > 1. For 
r 2 + z 2 <C 1, we require * — > (r 2 + z 2 ) -1 / 2 which is the star's intrinsic dipole 
field. Along the symmetry axis, \f(r = 0, z) = 0, and = 0) = 0. 

On the equatorial plane inside the light cylinder (r < 1, z = 0), B r = 
— (l/r)(d^/dz) = and H(r,z) = 0. For the closed field lines within the 
light cylinder ^ > ^ eq . On the equatorial plane outside the light cylinder 
(r > 1, z = 0), *(r,0) = * e(? = const and H(r,0) = 0. Thus for the open field 
lines, the range of is from zero on the axis to ^ eq on the equatorial plane. 

On the outer boundaries at r max and z max , we take free boundary conditions 
d 2 ^> /dn 2 = where n is the normal to the boundary. Other conditions have 
been tested including using equation (1) and they do not alter our results. 



2.2 Light Cylinder Condition 

A further condition on the solutions of equation (1) arises from the regularity 
of at the light cylinder. Because the coefficient of the A* term vanishes at 
the light cylinder, all field lines (\t values) which cross the light cylinder must 
have 

2^ = 4? = F for r = 1 and all z . (2) 
dr d^ 
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This relation determines H(^) for the field lines which cross the light-cylinder. 
In the following dH/d$l = H' . 
There are two possibilities: 

2.2.1 No Jets 

One is that all open field lines cross the light-cylinder This is the solution put 
forward by CKF. 

Because we calculate \& in a finite size region, J-(^>) is not determined for the 
open field lines which exit the region at z = z max inside the light-cylinder, r = 
to 1. These field lines have < * < * c , where * c = $(1, z max ). For the CFK 
solution, for this range of ^> we assume a linear interpolation, T = /<f C )F C , 
where the c-subscript indicates evaluation at (r = 1, z — z max ). The quantity 
T c is known owing to equation (2). Note that both ^ c and T c evolve as the 
iteration proceeds. 

For the CKF case, 

d^ T = -H 2 (^ eq -0+)>0 , (3) 

where + is an arbitrarily small positive quantity. Because H(*f> eq ) — 0, this 
requires a poloidal current sheet at z = + for r > 1 with 

J r = -^B r (r,Q+)H(y eq -0+)6(y -V eq + 0+) , 

= ~^(*^ - 0+)«(z - 0+) . (4) 

The poloidal current flow is sketched in Figure la. There is a corresponding 
positive current density oc S(z + + ) due to the lower half-space. Inside the light- 
cylinder the poloidal current sheet follows the dipolc-likc poloidal field lines to 
the star's surface. This is required in order to have H = within the closed field 
line region of the magnetosphere. At the same time there is a delta-function 
toroidal current flow with 

J - (c/47r)B r (r, 0+)S(z - 0+) (5) 

for r > 1. There is a corresponding positive toroidal current flow oc 6(z+0 + ) due 
to the lower half-space. Inside the light-cylinder there is also a delta-function 
toroidal current sheet associated with the mentioned poloidal current sheet. 

The magnetic field just above the equatorial plane for r > I is B + = 
B r i + B,p4> , and this is equal and opposite to the field just below the plane. 
The oppositely directed fields are unstable to annihilation. The electromag- 
netic stress T zz varies discontinuously from zero on the equatorial plane (by 
symmetry) to a finite value at z = + . 



L 
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2.2.2 Jets 



A second possibility is that there is a collimated jet along the +z axis (and — z 
axis) so that not all open field lines cross the light cylinder. The open field lines 
which do cross the light cylinder must obey equation (2). For the other open 
but collimated field lines (0 < r < tl, z — > oo), H(^f) is arbitrary. 

Thus H(^) is not determined for tp = */* c < 1. Because H(0) = we 
consider the simple dependence 



where kn is an adjustable constant. As kn increases in equation (6) the axial 
current flow within the light-cylinder (oc dH/dty) increases. Owing to equation 
(2), (HH')c is fixed at each step of the numerical iteration of 'J. Of course, 
both \I/ C and (HH') C evolve as the iteration proceeds. Consequently we can 
calculate 



at each iteration. 

We seek solutions without an equatorial current sheet. That is, we search 
for solutions with 



Figure lb shows a sketch of the poloidal current flow. 

3 Numerical Solutions 

The numerical calculations of $ for both cases of no jets and jets were done 
using successive iterations on a uniform 500 x 500 mesh in a square region 
r max x z max = 5x5. We have obtained the same results for 600 x 600 grid 
for a 6 x 6 region. The initial ^(r, z) used to start the iteration consists of the 
vacuum dipole field inside the light-cylinder and the straight line extrapolation 
of the field lines outside this cylinder. Convergence of the iterations is measured 
by the change of ^ between iterations. Figure 2 shows the solution ^(r, z) and 
H{^>) for the CKF case where there is no jet. 

3.1 Solutions with Jets 

Figure 3 shows \&(r, z) and H{^>) for the case of a jet flow along the axis and 
a disk-wind in the equatorial plane where kn — 2.48. The flow has zero net 
poloidal current flow, J** e9 (M>T = and no current sheets. The main param- 
eters of this solution are: along the z— axis "J = = H; on the light-cylinder 
(at z — 5) the values are *ff c and H c = H(^ c ); in the B p = region and 
Hoo] and on the equatorial plane ^ eq and H($f eq ) = 0. The numerical values 




(6) 





(7) 
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are * c = 0.0384, H c = 0.6180, = 0.133, = 0.637, and ^ eq = 0.318. 
This solution is not unique. For example, we have found an analogous solution 
for kn ~ I- 

The total power output into both the upper and lower halfspaces is E to t = 
B o r l c Jo'" dV[-H(V)] = 0.152So, where E = B 2 r 2 L c. In contrast, the to- 
tal power output of the quasi-spherical wind solution is E tot w B$r\c (e.g., 
McKinney 2006). 

The jet flow consists of a region of collimated flux within the light-cylinder 
r L where ^ < \I/ C = 0.0384, and a quasi-collimated region outside r L where 
* c < * < = 0.133. The power flow in the collimated jets is E c - ]ct = 

E /* c dtf [-#(#)] = 0.0155SO. The power flow in the quasi-collimated flows 
is £ qcjct = E J* 30 dtf [-.#(*)] = 0.0595£V The power flow in the disk wind is 
-Edwind = 0.0767£t)- Thus, about 10% of the total power goes into the collimated 
jet, 40% into the quasi-collimated flow, and 50% into the disk-wind. 



3.1.1 Radial Force Balance of Jet 



For conditions where a collimated jet exists, the z— derivatives in equation (1) 
vanish. The pulsar equation can then be written as 



i!.B 2 + — — 

dr z r 2 dr 



r 2 (Bl-E 2 ) 



= , 



(8) 



which expresses the radial force balance. Multiplying this equation by r 2 dr and 
integrating from the axis to r gives the radial virial equation 



H 2 = (r 2 



, dr /d*Y 

2 /„ " (* ) 



(9) 



following Lovelace, Berk, and Contopoulos (1991). 

Figure 4 shows the two sides of equation (9) obtained from our numerical 
solution at z = 5. There is a ~ 10% difference of the two sides of the equation. 
Calculations in a much larger region are required to determine whether or not 
the quasi-collimated flux outside the light-cylinder becomes collimated at very 
large distances. 



3.1.2 Vertical Force Balance of Disk 

The field solution shown in Figure 3 satisfies equation (1) accurately everywhere 
except close to the equatorial plane at r 2 3> 1. The reason for the discrepancy 
in this region can be understood by considering the vertical force balance for 
z 2 <C r 2 and r 2 » 1 where B 2 <C B 2 . For these conditions equation (1) is 
approximately 

2 a 2 * _ ldH 2 _ IdH 2 1 
~ r dl 2 " _ ~2~d¥ ~ ~2 dz d^/dz ' 
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This can be rewritten as 



= , (10) 

where £ = z/r. This expresses the vertical force balance near the equatorial 
plane: The term — (d^/dt;) 2 represents the negative pressure of the electric 
field and it gives an upward force. The H 2 term is magnetic pressure due to 
the toroidal magnetic field and it exerts a downward force. 

Figure 5 shows the vertical profiles of <J/ and H. These profiles do not 
satisfy equation (10). The reason is that equation (1) omits the plasma kinetic 
energy density in all space including z — where the magnetic field reverses 
direction. Therefore we include the influence of the kinetic energy density as 
a term K(^f)/2 within the square brackets of equation (10). The origin of the 
kinetic energy is from the annihilation of the magnetic field. [A term of this form 
can be derived from the relativistic Grad-Shafranov (GS) equation of Lovelace 
et al. (1986). The right-hand side of equation (86) of this work includes a term 
— Air pr 2 dJ {''Si) / d^> , where J = 7(1 — rv^/c) is the Bernoulli constant and 7 is 
the Lorcntz factor. We may write pr 2 = /(*) and dK/d^ = 8tt f(^)dJ ("$>)/ d*. 
The radial flow speed is v r = -cH^ 1 (d^ /d£) so that pr 2 v r = fct(*). ] Thus 
we obtain 

(H) -KiO + H'iO-Hi,, (11) 

where is the constant value of H at large £. Clearly it is necessary to have 
K + H 2 >H^. Because H(0) = and (<9*/<9£) 5=0 = 0, K{0) = H^. Because 
(d^/dt;) — > as £ increases, we have K — > as £. A sufficient condition for 
having E 2 < B 2 = B 2 + B 2 is K < H^. To illustrate the behavior we consider 
the dependences H 2 = H^[l - exp(-£ 2 /£| )] and K = cxp(-£ 2 /£ 2 r ) with 
S = Zk/Zh > 1. 

Figure 6 shows the £ dependences of \P, H, and rE z calculated from equation 
(11) for an illustrative case which maintains the conditions of Figure 5 of = 
\&eg — "J 00 = 0.185 and Hoo = 0.638. For this case we have taken 6 = 2. As 
a result, — 0.286 which corresponds to a half-angle thickness of the disk 
of ~ 16°. Near the equatorial plane the magnetic field has the form of an 
Archimedes' spiral, namely, dr/d(f> = — r^(l — (5 -2 ) 1 / 2 . With this modification 
of the disk configuration, the global field solution involves no delta-function 
current sheets. 

4 Conclusions 

This work describes a new solution for the structure of the magnetosphere of 
an aligned rotator described by the force-free pulsar equation. This is obtained 
by adjusting the current flow along the poloidal field lines which remain within 
the light-cylinder. When this current flow is sufficiently large a collimated jet 
forms within the light-cylinder and a quasi-collimated flow forms outside of it. 



d_ 

si 



1 /a* 

2 \d£ 



-H 2 
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The solution is not unique. The jet is collimated by the toroidal magnetic 
field. At the same time an anti-collimated disk-wind forms in the vicinity of the 
equatorial plane. The anti-collimation is due to the toroidal magnetic field. The 
vertical force balance of the disk-wind requires the inclusion of a finite kinetic 
energy density near the equatorial plane. Roughly one-half of the open field 
lines go into the jets and the other half to the disk-wind. The total current flow 
in the jets is equal and opposite to the current flow in the disk. Thus there is 
no current sheet inside or outside of the light-cylinder. 

Our jet/disk-wind solution represents an alternative to the quasi-spherical 
wind solutions where a major part of the current flow is in a current sheet. Such 
a current sheet is unstable to magnetic field annihilation. Furthermore, the way 
in which the configuration is reached is expected to be important. Consider a 
possible simulation experiment where the configuration is reached by dynamical 
evolution from an initially non-rotating star. The presence, distribution, and 
density of the initial plasma is important in that it is essential for the current 
flow as the star is spun up to a final rate. The value of the density of the 
background plasma was found to have an crucial role in determining the inflation 
magnetic loops threading a differentially rotating disk in relativistic particle-in- 
cell simulations (Lovelace, Gandhi, & Romanova 2005). The transition from 
one type of equilibrium to another remains to be investigated. 

We thank A. Spitkovsky, M. Sulkanen, and A. Timokhin for valuable dis- 
cussions. This work was supported in part by NASA grant NAG 5-13220 and 
NSF grant AST-0307817. 
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Figure 1 : (a) Sketch of the poloidal current flow in a pulsar magnetosphere with 
a quasi-spherical wind with no jet as proposed by Contopoulos et al. (1999). 
(b) Sketch of the poloidal current flow with a collimated jet and a disk-wind. 
This solution is related to that proposes by Sulkanen and Lovelace (1990). 




Figure 2: Quasi-spherical wind solution similar to that of Contopoulos et al. 
(1999). The jump in H at z = is related the equatorial poloidal current 
current sheet shown in Figure la. 




Figure 3: Pulsar magnetosphere with jet obtained by successive iteration of 
equation (1) with the constraint that f^*" d^>T = 0. H(*$>) within the light 
cylinder is given by equation (6) with kn = 2.48 and (3 — 0.9995. 
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Figure 4: Two sides of equation (9) obtained from our numerical solution at 
z = 5. RHS denotes the right-hand-side of the equation. 
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Figure 5: Vertical profiles of "J" and H at r = 4. At large z, ^>oo = 0.133 and 
iJoo = 0.638. 
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Figure 6: Vertical profiles of 'J, H, and rE z calculated from equation (11) for 
illustrative conditions but maintaining A^P = ^f eq — = 0.185 and Hoo = 
0.638. 
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